Neutron stars with small radii - the role of A resonances 
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Recent neutron star observations suggest that the masses and radii of neutron stars may be smaller 
than previously considered, which would disfavor a purely nucleonic equation of state. In our model, 
we use a the flavor SU(3) sigma model that includes A-resonances and hyperons in the equation 
of state. We find that if the coupling of the A resonances to the vector mesons is slightly smaller 
than that of the nucleons, we can reproduce both the measured mass-radius relationship and the 
extrapolated equation of state. 
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I. INTRODUCTION 

The mass-radius relationship of neutron stars is of key 
importance if one wants to understand the high-density 
low-temperature region of the hadronic equation of state. 
Depending on this relationship, certain models for the 
hadronic equation of state can either be confirmed or 
ruled out. 

A recent analysis of the masses and radii of neutron 
stars by Ozel et al. [U-Il] seems to imply that they are 
smaller than previously assumed, with radii between 8 
and 10 kilometers and masses between 1.6 and 1.9 so- 
lar masses. This information is extracted from the ther- 
monuclear bursts of x-ray emitting neutron stars in bi- 
nary systems. Using the measured distance, the appar- 
ent surface area during the cooling phase of the bursts 
and the flux in certain types of bursts which exceed the 
Eddington limit, a constraint for the masses and radii 
of the measured neutron stars is deduced. From this, 
a backward-extrapolation of error bars for equations of 
state [BHH that are in accordance with the measured val- 
ues is performed. A comparison with various models of 
equations of state suggests that a purely nucleonic equa- 
tion of state is too stiff, i.e. produces too high masses at 
too large radii. It is therefore possible that more exotic 
particles are realized within the interior of neutron stars 
[0| , which soften the equation of state. 

Several attempts have already been made to under- 
stand the measured values. Fattoyev et. al [l(| use RMF 
models with a nucleonic EoS. While they can fit the EoS, 
they are not able to fully reach the small masses and radii 
reported by Ozel et al. . Ciarcelluti et al. [ll| explore 
the possibility of neutron stars having a dark matter core. 
Drago et al. [T3] propose that this could be an experi- 
mental signal for the existence of two types of stars, or- 
dinary neutron stars and "ultra-compact" neutron stars 
like quark stars. Steiner et al. [l3[ performed a reanalysis 
of the experimental data with slightly different assump- 



tions. They relax the assumption of the stellar radius at 
maximum temperature being equal to the photospheric 
radius, increasing the allowed radius interval. As a re- 
sult, they reproduce mass-radius relationships which are 
in agreement with a purely nucleonic equation of state. 

In our approach we want to explore the possibility of 
additional particle species being present, namely A res- 
onances and hyperons, and investigate the effect on the 
equation of state of neutron stars. 



II. DESCRIPTION OF THE MODEL 

A hadronic SU(3) sigma-omega model is used to de- 
scribe the properties of neutron star matter. The 
included baryonic degrees of freedom are the baryon 
octet (n,p, and hyperons A, E + ' 0,_ , S 0,_ ), the leptons 
(e, n), and the resonances from the spin 3/2 decuplet 
( A ++ ' + <°>-, £*+> --, and Q~). The mesons, that 

mediate the interactions between the baryons, are the 
vector-isoscalar u and </>, the vector-isovector p and the 
scalar-isoscalar a and £ (strange quark-antiquark state). 

The Lagrangian density reads 
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Lxin is the kinetic energy term for the hadrons and 
leptons. In addition there is an interaction term between 
the baryons and the scalar and vector mesons 
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Lint = — Z)i V'i[7o(sww + gi^cj) + gz P T 3 p) + M*]ipi 

with the effective mass M* given by 

M* = g la <j + g l( ( + 6mi, (3) 
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with a small bare mass term Srrii. The coupling strengths 
of the baryons to the scalar fields are connected via SU(3) 
symmetry relations and the different SU(3) invariant cou- 
pling strengths are fitted to reproduce the baryon masses 
m vacuum (seeRef. [U| for the values of the couplings). 

The self-interaction terms for the scalar and vector 
mesons read, 



Lseif = + m 2 p p 2 + m 2 (j) 2 ) 



+fc ( ( 7 2 + C 2 ) + fc 1 (a 2 +C 2 ) 2 
+k 2 + C 4 ) + k 3 a 2 C + k 4 In (4) 

and the explicit chiral symmetry breaking term is given 
by 
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In the case of the baryonic octet we use an f-type 
coupling between the baryons and vector mesons, which 
yields coupling strengths as given by quark counting 
rules, i.e. g iLU = (n* - n\)g% , g lcf , = -{n\ - n\)y/2g% , 
where g% denotes the vector coupling of the baryon octet 
and n l the number of constituent quarks of species i in a 
given hadron. A more detailed description of the model 
can be found in [T3|. The model has been tested ex- 
tensively both for nuclear matter [TBI and neutron stars 

The parameters are fitted to reproduce the correct vac- 
uum masses of the baryons and mesons, nuclear satura- 
tion properties at po = 0.15 particles per fm 3 (the cor- 
rect binding energy per nucleon, B/A = - 16 MeV and a 
compressibility of K = 297.32 MeV) and the asymmetry 
energy (E sym = 32.5MeV). 

As the values of the vector coupling of the decuplet 
are not constrained by the properties of saturated nuclear 
matter as in the case of the baryon octet, in the following 
we allow for moderate deviations from the overall vector 
coupling strength for the baryon decuplet compared to 
the octet. In order to study the effects of such a deviation 
in a most direct way we introduce a single parameter r v , 
which in the case of the A resonances is defined as 
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9Nu 
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with the same rescaling of the other states of the de- 
cuplet. 

Using this approach we solve the equations of mo- 
tion by differentiating the grand-canonical potential S = 
— p = e — pqb of the model for T = with respect to the 
fields (in mean-field approximation) . Then we determine 
the equation of state and use it to solve the Tolman- 
Oppenheimer-Volkoff equations for spherical static stars. 
A standard BPS crust EoS |l8| is added to the outer 
layers of the star. 
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FIG. 1: Mass- Radius-relationship of Neutron stars for various 
couplings of the A resonances, starting from r v = 1 (upper 
line) to 0.8 (lowest line). Also included are the 1-a errorbars 
for measured neutron stars from [l|. The black diamond on 
each curve represents the maximum stable configuration of 
the neutron star. 



III. RESULTS 

A resonances have a repulsive vector potential which 
works to counteract gravity in a compact star. If the 
interaction strength of the A resonances with respect to 
the nucleons is lowered, i.e. r v < 1, then the equation of 
state becomes softer and hence the maximum supported 
mass of the star decreases. 

Figure 1 shows the impact of various coupling strengths 
r v of the A resonances on the overall mass-radius rela- 
tion of the neutron star. The upper black curve repre- 
sents r v = 1, and in each subsequent curve the coupling 
strength is lowered by 0.05, arriving at 0.8 in the lowest 
curve. If the coupling strength of the A resonances is re- 
duced, the maximum mass of the star likewise decreases. 
We found that for values of 0.7 or lower, the equation 
of state has to be considered with care and may start to 
yield unphysical results such as the appearance of As in 
the nuclear ground state. 

However, for values of r v around 0.9, which corre- 
sponds to the physical case where the coupling strength 
of the A resonances is only slightly lowered compared 
to normal nucleons, the mass-radius-relation is in good 
agreement with the values given by Ozel et al. . 

Figure 2 shows the equation of state for various cases 
considered. For comparison reasons, we use the same 
units as Ozel et al. [l| and include their error bars (black 
perpendicular lines) . The stiffest line represents a purely 
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FIG. 2: Considered equations of state(from stiffest to softest): 
purely nucleonic (solid), including A resonances at r v = 1 
(dashed), A resonances and hyperons at r v — 1 (dotted), A 
resonances at r v = 0.9 (dashed), A resonances and hyperons 
at r v = 0.9 (dotted). The additional inclusion of hyperons 
only creates a minor effect. 
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FIG. 3: Particle abundancy as a function of the baryonic den- 
sity in l//m 3 for the case with hyperons and Delta resonances 
present, rv =0.9 . 



nucleonic equation of state, n, p and e being present. In 
all subsequent cases A resonances are included. We look 
at two cases, r v = 1 and r v = 0.9 (both dashed). The 
dotted lines below each of the curves represents the effects 
of muons and hyperons being added. The inclusion of ad- 
ditional particle species generally softens the equation of 
state. The major effect comes from the delta resonances, 
while the hyperons and muons only make minor changes. 

Figure 3 shows the particle abundances for the best-fit 
case of r v = 0.9, with nucleons, leptons, A resonances 
and hyperons present. At densities between 0.2 and 0.5 
particles per fm 3 , where go = 0.16/m' — 3) represents 
normal nuclear matter density, the A resonances start to 
appear, first the A - then the others. At very high densi- 
ties beyond six times qq, the hyperons start to appear. In 
principle all the listed particles are present in the model 
and can appear, including the excited spin 3/2 states. 
However, as one can see in Figure 3, only the A starts to 
appear in a small fraction in the relevant regime. 

It is noteworthy that negatively charged particles like 
the muon and the A~ at a certain point take up the role 
of the electrons in charge conservation and thus reduce 
the electron abundance. At a certain point, the A par- 
ticles, due to their reduced repulsive potential, become 
the dominant particle species. 

Figure 4 shows the effective mass of the nucleons, A 
and A particles within the chiral model as a function of 
the density. 
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FIG. 4: Effective mass of the particles as a function of the 
baryonic density. 



IV. DISCUSSION 

Recent experimental data [l[ seem to favor neutron 
star masses and radii that are smaller than previously 
considered. An equation of state consisting only of nucle- 
ons and electrons is not able to reproduce the measured 
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masses and radii. 

Note, however, that in a different analysis, Steiner et 
al. [l3[ relax the assumption of the stellar radius at maxi- 
mum temperature being equal to the photospheric radius, 
increasing the allowed radius interval. As a result, they 
reproduce mass-radius relationships which are in agree- 
ment with a purely nucleonic equation of state. If the 
analysis of Oezel et al. is correct, and neutron stars 
have smaller masses and radii than previously consid- 
ered, then this may hint to more exotic particle species 
being present in the equation of state. We have shown 
that an equation of state with A resonances in such a 
case would be in accordance with the experimental data. 
While our approach includes hyperonic degrees of free- 
dom in a natural way, the results are dominated by the 
impact of the A resonances, whereas hyperons only have 
a minor effect on the result. 

If, on the other hand, Steiner et al. should be correct 
with their analysis, this would point to a more conser- 
vative mass radius relation and equation of state, which 
can be described by the presence of nucleons and elec- 
trons alone. A greater amount of analyzed experimental 
data may be necessary to resolve this issue. 



V. CONCLUSIONS 

We have applied a non-linear realization of an extended 
sigma-omega model, including A resonances and hyper- 
ons in the equation of state of a neutron star calcula- 
tion. A parameter scan for the coupling strength of the 
A resonances has been performed. By slightly reducing 
the strength of the vector coupling by about ten percent 
compared to the values for the baryonic octet we obtain 
good agreement with the results from Ref. both in 
the mass-radius relation and equation of state. It is note- 
worthy that this agreement can be achieved by the not 
too-exotic inclusion of A resonances in the equation of 
state without considering quark cores or quark stars. 
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